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Nomenclature 


A = upstrean location of pressure tap 
kj = flow normal area = V^,/L 
A„ = speclffieii normal area 

= wetted wail area (total wail area exposed to fluid) 

B = downstream location of pressure tap 
= specific heat at constant pressure 
Dj, - specimen hydraulic diameter 

f = friction factor 
f = normalized heat flux on target 
fq(x,y,v) = heat flux distribution function 
G = mass flow rate per unit flow normal area 
= m/Aj = pV 

h = heat transfer coefficient 
h = enthalpy 

j(x,y,v) = normalized heat flow through heat flow meter 
= heat flow through heat flow meter 
k = thermal conductivity 
K = heat flow meter calibration constant 
L = length of specimen 
m = mass flow rate 
Nu = Nusselt number = h-D^/k 
P = pressure 

Pr = Prandtl number = /i-Cp/k 

q(x,y,v) = local heat flux based on area normal to the furnace target 
= heat leak to manifolds through insulation 
^'in " beat leak to inlet manifold 

Qp, - fraction of total heat flow on specimen added up to position x 
= integration of furnace calibration function f ^ , 0 to x 
Qy - total heat transfer to spec^ n 

- local heat flux (heat flow per unit area) into the cooling fluid based on 
total wetted-wall area of the specimen 
r = recovery factor = for turbulent flow 

Re = Reynolds number = pV\/^ 

T = temperature 

” cooling fluid adiabatic wall temperature 
Tf = local bulk fluid temperature 
= Eckert reference temperature 
= specimen wall temperature 

V = heater voltage 

V = velocity 

V(x,y,v) = measured heat flow meter vol^’age 
Vo = open volume in specimen 
w = heat flow meter half-width 

X = position coordinate parallel to flow direction 
Xg furnace target half -width » L/2 

y » position coordinate perpendicular to flow direction 

P s coefficient of thermal expansion 
o = standard deviation 
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ft = dynamic viscosity 
p = density 



0 * location In inlet maulfold of PRT 

1 = location in outlet manifold of PRT 


vlli 


i 



Apparatus for Measuring High- Flux 
Heat Transfer in Radiatively Heated Compact Exchangers 


Douglas A. Olson 

Chemical Engineering Science Division 
National Institute of Standards and Technology 
Boulder, CO 80303-3328 


We describe an apparatus which can deliver uniform heat flux 
densities of up to 80 W/cm^ over an area 7.8 cm by 15.2 cm for use 
in measuring the heat transfer and pressure drop in thin (6 mm or 
less), compact heat exchangers. Helium gas at flow rates of 0 to 
40 kg;Ti and pressures to 6.9 MPa (1000 psi) is the working fluid. 

We describe the instrumentation used in the apparatus and the methods 
for analyzing the data. We will use the apparatus initially to 
test the performance of prototype cooling jackets for the engine 
struts of the National Aerospace Plane (NASP) . 

Key words: apparatus; compact heat exchanger; convection heat 
transfer; friction factor; high temperature; National Aerospace 
Plane; pln~fin; radiative furnace. 
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1 Introduction 


Research at NAS on vehicles for hypersonic flight has been active since 
the 1960* s (Shore, 1986). The renewed interest in developing a National 
Aerospace Plane (NASP) requires solving many demanding heat transfer problenis. 
associated with the anticipated hlgii heat fluxes on va: us components. The 

problem of interest motivating this work is the cooling of the engine struts 
Due to aerodynamic heating associated w^th the combustion of the hydrogen 
fuel, along with radiation from the fuel comDustion, the engine struts are 
expected to receive a heating load in excess of 200C W/cm^ (Scotti et al., 
1988). NASA plans to cool the struts by attaching a cooling jacket heat 
exchanger to the surface facing the high heat flux. Hydrogen gas will flow 
through the cooling jacket, absorbing the heat before entering the engine. 

The anticipated conditions are Chat the hydrogen gas will enter the heat 
excha^^gers at 56 K and 6.9 MPa (1000 psi), and exit at 890 K and 9.8 MPa (700 
;ri). The heat exchangers are expected to be thin perpendicular to the flow 
direction (6 nun or less) to add minima] weight and t. ickness tj the struts. 
Small flow passages will also produce high rates of convective heac transfer, 
which will reduce the exchanger temperatures. Reynolds numbers are expected 
to be in the range 10000 to 30000, with the variation due to the flow rate 
and the specific design of the flow passage. Because of the small flow 
passages and high heating rates, the heat transfer performance of the 
prototype heat exchangers is beyond the range of past experimental work and 
correlations . 

A leading candidate for the cooling jacket is the **pin-fin'' heat 
exchanger (fig. 1). This design is being considered based on fabrication 
techniques, structural considerations , and the anticipated heat transfer 
performance (Buchmann et al., 1985). Nickel 201 is one of the pr posed 
materials, due to its high ductility and high thermal conductivity. A pin-fin 
heat exchanger consists of a closely spaced channel with pins connecting the 
top aiiJ bottom surfaces. The pin-fins increase the heat transfer area; by 
promoting .^rbulence they also Increase the local heat transfer coefficient. 
Compared to other compact heat exchangers, the pin- fin configuration has high 
heat transfer per unit volume, but also a large pressure drop (Kays and 
Lx>ndon, 1969). Pin- fin heat exchangers have been used for cooling turbine 
blades in gas turbines, and most of the existing correlations are based on 
work for chat application (Armstrong and Winstanley. 1987). 

The objective of this work was to design and '?uild an experimental 
apparatus which could test prototype heat exchangers at some of the conditions 
anticipated for the NASP application. Test results would be heat transfer 
performance and pressure drop as a function of Reynolds number, heating rate, 
and flow geometry. The apparatus will be used initially to test pln-fln heat 
exchangers, and In the future will be capable of testirg any flat, compact 
heat exchanger which can conform to the furnace dimensions. 

2, Description of the experimental apparatus 

The experimental apparatus has been designed to initially test a subset 
of the conditions required for the NASP application. Those conditions are (1) 
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a healing rate of 0-80 W/crn^ ; (2) an Inlet temperature c f ^00 K H) a 
cooling-gas pressure of 6.9 MPa at the Inlet, and (4) an outlet temperature of 
810 K. The size of specimen which can be tested is 7 8 err vide bv 15 2 rm 
long Because of the stran^ritics in specific neat, Thermal conduc r 1 vi t v , and 
dynamic viscositv between helium and hydrogen, helium will be the first 
cooling gas tested. The Reynolds numbers and temperature rise from specimen 
inlet to outlet can be matched between helium and hydrogen As the taking of 
data progresses, indicating the critical experimental parameters, flow loop 
components can be added which yield cryogenic inlet conditions or higher 
heating rates. Hydrogen gas will be tested after moving the apparatus to a 
special site for hazardous tests and adding hydrogen- safe ty features. 

2.1. Flow loop 

The helium flow loop is shown in figure 2, with the details of the 
specimen furnace section in figure 3. Helium gas at 13.8 MPa (2000 psi; is 
supplied from a tube trailer outside the laboratory. T^e tube trailer 
contains 1100 m^ of gas (STP). Valve #1 at the trailer isolates the gas from 
the 1 . 2 cm diameter high pressure copper piping to the laboratory A second 
shutoff valve (#2) within the laboratory isolates the inlet piping from the 
portion within the laboratory. With valves «1 and «2 open, gas flows from 
trailer, through the inlet piping, and is filtered before entering ^he dome 
loaded pressure regulator (valve #3). The regulator sets the flow pressure 
downstream of the regulator to the value of an external control pressure. 6 ^ 
iir« viuOO psi) or less. A pressure relief valve downstream of the regulator 
set at 8.3 MPa (1200 psi), prevents system over pressure The gas flows 
through stainless steel tubing. 1 9 cm OD and 0 21 cm wall, from the pressure 
regulator to the furnace . 

Within the furnace (fig. j), the gas flows into an inlet di s t r ibut i ori 
manifold (347 stainless steel) which directs the gas to the heat exchanger 
specimen. A similar distribution manifold collects the ^as exiting the 
specimen and directs it to the outlet tubing. The heat exchanger specimen is 
brazed into the space between the manifolds. The specimens are loc::trd in the 
target area of the furnace, which deli radiant heat to the specimens and 

raises the temperature of the helium as it flows through the specimen 

The furnace consists of a high- intensity radiant infrared heater, 
surrounded by refractory insulation, 5 cm thick or greater, which directs and 
re-radiates the heat from the heater to the specimen. The heater contains six 
high' temperature infrared lamps mounted in an aluminum housing Filtered shop 
air at 0.8i m^/min (30 cfa) or greater cools the heater lamps, and water at 4 
kg/m In (I gpm) ot cools the housing. 77ie radisnt heater is powered by 

a phase-angle power controller which uses 480 VAC. .single phase, and 75 A ar 
full pow»r setting. At the naxiaua heat flux setting, the refractory walls 
reach an aatiaated teaperature of 1760*C n200'F). 

Downatraaa of the furnace section, the hot gas flows through aore 
stainless steel tubing (1.9 ca OD, 0.21 ca wain to a cooling coll (, 5 a or 
the saac tubing) ianersed in a Katar bath. Tf* ' water bath cools the he.*. urn 
temperature to 30*C or less. The rate of gas flow is manually adju;jted at 
the bath outlet by valve a4, which also drops the gas pressure to atmospheric 
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pressure. Beyond the valve, helium flow rate is measured with a heated- tube 
thermal mass flow meter, with ±1% uncertainty in the reading (at a 95% 
confidence level) . After exiting the flow meter the gas flows into an exhaust 
line of 5 cm copper tubing, and is vented outside the laboratory. 

2.2. Instrumentation 

Measurements are made of temperature of the gas in the inlet and outlet 
manifolds, gas pressure in the specimen, specimeti temperatures, and the 
aforementioned gas flow rate. Uncertainties in the measured quantities are 
summarized in Table 1. 


Table 1. Uncertainties in experimental measurements and gas properties 

at a 95% confidence interval 


Measurement /Proper tv 

Technique 

Major Source 
of Uncertainty 

Magnitude 
of Uncertainty 

Gas Flow Rate 

Thermal Mass Flow 
Meter 

Meter Calibration 

±1% 

Heat Flux 

Calibration cf 
Furnace 

Heat Flow Meter 

±4% 

Gas Inlet and 
Outlet Temperatures 

Platinum Resis- 
tance Thermometer 

Radiation 

±0.5 K 

Gas Pressure 

Pressure Trans- 
ducer 

Calibration 

±0.25% of 
full scale 

Specimen Temperature 

Type-N Thermo- 
couple 

Wire Calibration, 
Installation 

greater of 
±0.4% or 1.1 K 

Gas Temperature in 
Specimen 

Integrate Energy 
Equation 

Heat Flux 
Uncertainty 

±1% of Tj-To 

Gas Density 

Thermodynamic 

Function 

Function Accuracy 

±0.1% 

Gas Enthalpy 

Thermodynamic 

Function 

Function Accuracy 

±0.2% 

Gas Specific Heat 

Thermodynamic 

Function 

Function Accuracy 

±5% 

Gas Viscosity 

Thermodynamic 

Function 

Function Accuracy 

±10% 

Gas Thermal 
Conductivity 

Thermodynamic 

Function 

Function Accuracy 

±3% 
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We determined the heat flux distribution on the specimen by calibrating 
the furnace prior to inserting the specimen. The heat fliix distribution is 
defined as the local heat flux as a function of position over the furnace 
target and as a function of heater lamp voltage. Details of the furnace 
calibration are given in appendix A. The uncertainty in the heat flvix is 
estimated as ±4%. 

The gas inlet and outlet temperatures are measured with platinum 
resistance thermometers (PRTs), 0.48 cm diameter, Inserted in the gas 
manifolds at locations 0 and 1 of figure 3. The uncertainty In the 
thermometer calibration is ±0.08 K; the uncertainty introduced by the digital 
voltmeter is ±0.036 K. Radiation errors in the thermometers are estimated as 
less than ±0.5 K, while conduction errors in the thermometers are estimated as 
less than ±0.01 K. The gas pressure at the specimen inlet (location A in 
figure 3) is measured with a variable -reluctance pressure transducer. 

Pressure difference between the specimen inlet and outlet (A to B) is measured 
with a second pressure transducer, also a variable-reluctance tjrpe with 
interchangeable diaphragms. The uncertainty in the calibration of both is 
±0.25% of full scale. Uncertainty introduced by the voltmeter is ±0.011% of 
the pressure reading. 

Specimen temperatures will be measured with type-N thermocouple wire, 
wire diameter 0.25 mm. An array of thermocouples, 25 or more depending on the 
specimen design, will be spot-welded to the side opposite the radiant heat 
flux (insulated-side). Temperatures measured at the insulated-side 
thermocouples will be used to determine the heat transfer coefficient. The 
relative uncertainty from pair-to-pair is ±0.4% of the reading or ±1.1 K, 
whichever is greater. At the maximum expected specimen temperature this is 
±2.1 K. Estimates of the heated-side temperature will be made with type-N 
thermocouples mounted in a fashion shown in figure 4, which is the method for 
the pin-fin specimen. Each wire of the pair is spot-welded to the heated 
surface, with the lead extending out a hole drilled through the specimen. For 
the pin-fin specimen, the holes will be drilled through the center of a pin. 
The theirmocouple circuit is completed by the specimen wall between the two 
wires. A quartz sleeve, 0.5 mm diameter, electrically Insulates the wire from 
the wall of the hole. Because a portion of the pin is removed and replaced 
by wire plus quartz, each of which has a lower thermal conductivity than the 
pin (nickel) , a pin with a wire mounted in it will be hotter than a pin with 
no wire mounted in it. We estimated the magnitude of this temperature rise 
from a finite-element analysis as 2-5 K at a radiant heat flux of 50 W/cm^ . 

All type-N thermocouples will be connected to an isothermal reference 
box. We measured the temperature uniformity of the reference box as ±0.05 K 
or better. Copper conductor wire runs from the reference box to the data 
scanner. The temperature of the reference box is measured with a platinum 
resistance thermometer (PRT) . Using the box PRT sensor for the reference 
temperature, we compared temperatures indicated by thermocouples connected to 
six terminals in the reference box with temperatures Indicated by the gas 
inlet and outlet PRTs. All sensors were immersed in a stirred water bath at 
20-30*C. This provided a system accuracy for the reference box pl\js PRT. 

The agreement in temperature between the thermocouples and gas PRTs was better 
than ±0.06 K. 
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All instrument signals are multiplexed through an automated scanner and 
measured on a digital voltmeter. The scanner and voltmeter are controlled 
with a personal computer through an IEEE 488 bus. Raw signals are stored on a 
hard disk and copied to floppy disk for backup. Signals are converted to 
engineering units and the data analyzed at the completion of an experimental 
run. Some signal readings are converted immediately to engineering units and 
displayed on the video terminal to assist in monitoring and operating the 
experiment. The measurement uncertainties introduced by the data acquisition 
system have been included in the stated uncertainties of each sensor. 

3. Methods of data analysis 

The data required to evaluate the performance of the heat transfer 
specimens are the heat transfer coefficient, h, and the friction factor, f. 

3.1. Friction factor 

The friction factor results from an integration of the one-dimensional 
momentum equation in the flow direction: 

P^-Pb = + (2GVDh)£(f/p)dx, (1) 

where P - pressure; 

G - mass flow rate per unit flow normal area 

= m/Af pV; 

Af = flow normal area = open voluire/specimen length; 
p = density; 

V = velocity; 

= specimen hydraulic diameter; 

A = upstream location of pressure tap; 

B = downstream location of pressure tap. 

The first term on the right hand side of the equation is the pressure change 
due to flow acceleration, and the secona term is the pressure drop due to 
frictional effects. For the pin-fin specimens, the hydraulic diameter is 
given by (VanFossen, 1982): 

Dh - 4V,/A,, (2) 

where = open volume in specimen; 

A^ = wetted wail area (total wall area exposed to fluid). 

This definitior is ^appropriate when the pin height, pin diameter, and pin-tc- 
pin separation are of the same order, as all are likely significant In 
generating turbulence which drives the pressure drop and heat transfer. 

If the density change is small compared to the absolute density, and the 
pressure drop through the specimen is llne^-r, then the Integral can be 
approximated as a constant and the resulting equation for f is: 
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( 3 ) 


f = Pa-Pb - G2(1//>,-1/x>a) 

2(GV/>)(VDh) * 

with p = (/>A + Pa)/2* 

The density change criterion is met when the gas flows through the 
specimen without being heated, but when the specimen is heated the exit 
density can be less than half the entrance density and eq (3) is not valid. 
The friction factor is calculated for a series of tests with flow rate as the 
variable but with no heating, and a functional dependence of f on Reynolds 
number is developed, where the Reynolds number is defined as 

Re = pVDh/p. (4) 

Then, for experiments with specimen heating, the density and Reynolds number 
(evaluated using the local gas temperature and pressure) are calculated as a 
function of x. The integral in eq (1) is evaluated to calculate the pressure 
drop from A to B, which is checked against the measured pressure drop. This 
provides a technique for evaluating If f depends on more than Reynolds number 
alone. 

3.2. Heat transfer coefficient 

The heat transfer coefficient, h, is defined through the equation 

q. = h.(T,-T,„), (5) 

where q^ = local heat flux (heat flow per unit area) into the cooling 
fluid based on total wetted-wall area of the specimen; 
h = heat transfer coefficient; 

T^ = specimen wall temperature; 

T*w = adiabatic wall temperature of the cooling fluid. 

The adiabatic wall temperature is used in gas flows whenever the kinetic 
energy is significant compared to enthalpy changes (Rohsenow and Choi, 1961). 
Friction can cause the local wall temperature to exceed the bulk fluid 
temperature for an adiabatic specimen, and the adiabatic wall temperature 
approximates this effect. It is defined as 

T.w = Tf + rVV(2Cp), (6) 

where T^ = local bulk fluid temperature; 

r =5 recovery factor = Pr^'^ for turbulent flow. 

Adiabatic heating can be on the order of 1 K for the anticipated conditions. 
The local heat flux in eq (5) is expressed in terms of the total heat 
transfer to the specimen, , the total wetted wall area, and the furnace 
calibration function, f^, which is a dimensionless expression of the local 
normal heat flux: 
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The wall temperature used in eq (5) is measured with a thermocouple on the 
insulated side of tb^ specimen. 

In the defining equation of h we assume (1) that the heat flux is uniform 
into the cooling fluid at a particular location (the same at the top and 
bottom surfaces and at the pins Joining the top and bottom), (2) that the wall 
temperature is uniform, and hence (3) that the heat transfer coefficient is 
constant. However, since the specimen is heated from one side only and heat 
must flow through the specimen structure ^o enter the fluid from the pin and 
insulated wall, wall temperatures will not be constant. The dimensions of the 
specimen prevent instrumenting around the perimeter of a pin, and imbedded 
thermocouples alter the specimen temperature distribution, so the insulated- 
side thermocouple is used as the wall temperature. A finite-element 
conduction analysis using anticipated values of the heat transfer coefficient 
and the pin-fin configuration of figure 1 indicated that the temperature along 
the insulated side varies by 0.1 K or less with respect to pin position (for 
50 W/cra^ hot side heat flux). In addition, the analysis indicated that the 
temperature of the insulated side is very close to the fluid-wall interface 
temperature along the bottom surface (within ±1 K for a pin-fin specimen with 
50 W/cm^ hot side heat flux) . This compares with a temperature difference 
between the wall and bulk fluid of at least 10 K for the same conditions. 

These conclusions could change for specimen geometries and materials different 
from those of the pin-fin configuration. 

Combining eqs (5), (6) and (7) and rearranging, we get 

h . (8, 

(T.-(T[ + (rV“)/(2c7)ll 

We will calculate the heat transfer coefficient at all locations where 
temperature of the insulated side is measured with a thermocouple. We will 
measure the wall temperature directly, while the furnace calibration function 
fq is a function of location and heater voltage. The flow direction energy 
equation is used to calculate (to follow). Gas temperature T^ is 
calculated using the flow direction energy equation along with the furnace 
calibration (rlso to follow). 

The total heat absorbed by the heat transfer specimen equals the total 
heat absorbed by the specimen plus manifolds, less the heat leak through the 
furnace insulation into the manifolds Q;. . It is calculated from the 
temperatures of the gas Inlet and outlet, the gas pressure drop, and an 
estimation of the manifold heat leak. Or 

Qt = m vhj -ho) - Q.. (9) 

where h = enthalpy ; 

0 = location in inlet manifold of PRT; 

1 - location in outlet manifold of PRT; 

0^ - heat leak to manifolds through insulation. 

Heat absorbed by the manifold is typically 2-5% of the total heat flow, which 
we measured during calibrations on the furnace prior to inserting the test 
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specimens. We neglect kinetic energy changes from 0 to 1 since they are 
Insignificant compared to the uncertainties of the temperature measurement. 
The enthalpy change is given by f 

hi-h« - Cp (Ti-To) + £l(l-^)/p]dP. \(10) 

where = coefficient of thermal expansion. 


The pressure term is significant since helium is not an ideal gas at these 
temperatures and pressures. The integral is evaluated using the virial 
equation of state for the gas (McCarty, 1973) and is Included as a subroutine 
in the data analysis program. For that evaluation, pressure in the inlet and 
outlet manifold is required, although it is measured only at locations A and 
B (fig. 3), Since the pressure term is usually less than It of the 
temperature term, pressure is assumed to vary linearly along the specimen and 
to be constant within the end manifolds. The ratio of the specimen length 
between taps A and B to the length between A and 0 (the specimen inlet) is S. 
Thus, Pg and P^ are estimated as 

Po * Pa + <Pa-Pb)/8. 

Pi = Pb - (Pa-Pb)/8- (H) 

Combining eqs (9) and (10) yield for ; 

Qt * mlCp (Ti-Tg) + jj(l-fiT)/p]dP) - Q. . (12) 


The fluid temperature, T^ , is calculated by integrating the flow energy 
equation from the inlet nuinifold up to the locution of interest (designated as 
x), now Including kinetic energy: 


T„ * Tg 


Qi:.Sa. 

mc» 




£ 




}L1 

2c„ 


(13) 


where Qp, « fraction of total heat flow on specimen added up to 
position x; 

= integration of furnace calibration function f ^ , 0 to x; 
s heat leak to inlet manifold. 

The fluid temperature requires the velocity at x, given by 

V, =m/(A,p,), (14) 

and Che density is given by the equation of state (McCarty, 1973): 

» P«(T„,P,). (15) 

The pressure at x is estimated by 


9 



(16) 


P, » Pa - (PA-Pi) xA. 

With «q (16) sub«tltut«d into «q (13) to ovaliuite th« prevsur* tera, oqs 
(13). (14), and (15) fora • syatoa of throe equotlono in the unknowns of 
Ceapereture, velocity, end density. They ere solved throuj^ iteration by 
first neglecting velocity in eq (13) to calculete teapereture; eq (15) is used 
to celculete density, eq (14) is used to celculate velocity, end teapereture 
is re-calculeted froa eq (13) with velocity included. The iteretion continues 
until convergence (<0.01 K change in auccessive coaputetlons of teapereture). 
Only two or three iterations ere usually required since the kinetic energy is 
about 1% of the enthalpy change. 

With T, and V deterained at location x. the heat transfer coefficient is 
calculated using eq (I). The Nusselt nuaber, Prandtl nuaber, and Reynolds 
nuaber are then calculated: 

Ru * h'Dj|A> 

Pr » pCy/k. (17) 

The Eckert reference teapereture is used to calculate the fluid properties due 
to the large wall to fluid teapereture difference, as done by VanFossen 
(1982): 


T, » 0.5T, + 0.28T, ♦ 0.22T... (18) 

Fluid properties are calculated froa the functions given in McCarty (1973). 

4. Susoary 

Ve have described an apparatus which can aeasure the heat transfer and 
friction factor in thin, coapact heat exchangers which are radiatively heated 
in a furnace at heat fluxes of up to 80 V/ca^. Instruaentation and aethods of 
data analysis are described. The apparatus will be used initially to test a 
coapact heat exchanger of the pin-fin configuration using hellua gas as the 
heat transfer fluid. 


Ve acknowledge the assistance of Robert R. Wilson in perforaing the 
furnace calibration, developing date acquisition prograas, and constructing 
the apparatus. We appreciate the assistance of David E. Daney in designing 
the apparatus. The theraocouple aounting technique for the pin-fin 
configuration was developed by NASA Langley Research Center. 
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Appendix A, Calibration of the heat-flux distribution within the furnace 

Evaluation of the heat transfer performance of the specimens requires 
values for gas> temperature inside the flow passages and for the local heat 
flux incident on the specimen. Neither of these will be measured directly 
during a test run with the specimen in place; the local gas temperature can be 
inferred from an energy balance on the gas if the heat flux to the specimen is 
known. Variations in heat flux across tne fuma^t target area, the "heat flux 
distribution" are measured prior to the tests on the pin- fin specimens th*-ough 
tests on a "calibration specimen." For specimens with the same emlssivlty 
and thermal conductance to the heat sink, the furnace performance will not 
depend on the type of specimen receiving the heat flux. 

A.l. Description of the calibration apparatus 

To calibrate the furnace, the calibration specimen was placed in the 
target area of the furnace (see fig. A.l). Three heat flow meters soldered to 
the specimen measured the heat flow over the meter area, approximately 2.52 
cffl^ . With the radiant heater at a steady voltage and with water cooling the 
specimen, the plate with meters attached was traversed over the furnace target 
area to measure variations in Che x and y directions at that voltage. The 
readings from the meters were used only up to the point where the outer edge 
aligned with the furnace wall. The voltage was Chen changed and the x-y 
traverse was repeated. The measured meter heat flow was Chen analyzed to give 
a heat flux , and a function for the heat flux was developed with inputs x, y, 
and heater voltage. 

The calibration specimen consisted of a lower brass plate, 2.5 cm thick, 
and an upper copper plate, 0.64 cm chick. Cooling channels for water flow 
were machined in the brass piece. The plates were 12.7 cm wide and 33 cm 
long. The three heat flux meters were each 1.59 cm by 1.59 cm and 0.16 cm 
thick, and were located as shown in figure A. 2. The signal lead wires lay in 
shallow channels machined in the top copper plate, and copper shading shims, 
also 0.16 cm thick, covered the leads to prevent overheating. The top surface 
of the assembled piece, including shims and heat flow meters, was spray- 
painted a flat black to provide a uniformly high absorptivity. 

The calibration specimen was placed in r traversing mechanism which 
moved the specimen In the x and y directions in the plane of furnace target. 

In the x direction (parallel to the axis of the heater), the traversing 
mechanism could move all heat flow maters from one end of the furnace to the 
other. In the y direction (perpendicular to the axis of the furnace), motion 
waa limited to about 3.8 cm, sufficient to allow meters 2 and 3 to align with 
the +y wall, and meter 1 to align with the -y wall. This allowed about 1.9 cm 
rf meter overlap in tha center; however each meter did not cravetse the full y 
width. Position was read to 10.5 mm. At maximum specimen displacement in the 
traversing mechanism (meters aligned with the furnace wall), coverage of the 
target by the specimen was lOOt. 

The temperature of the calibration specimen was monitored with a type K 
thermocouple placed under the left-side shim of figure A. 2. Type X 
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charnocouples itonltored teaperaturas at an Intarior location of tha furnaca 
Insulation, on tha housing of tha haatar, and in the cooling air as It exited 
tha heater. The haatar voltage was also neasurad (as a percent of autxlmun 
voltage) . The cooling water flow rata was nonltored with a rotaaietar 
flowmeter, and the water temperature rise across the specimen was measured 
with thermistors. 

A. 2. Operation of the heat flow meters 

The heat flow meters generated a DC voltage proportional to the heat flow 
through them. They consisted of three layers of metal; the top and bottom 
were the EP thermoelem''nn, thermoelectrical ly different from the EN 
thermoelement in the middle layer. Flow of heat through the meter generated a 
temperature difference across the middle layer. Due to the alternating 
materials of construction, the temperature difference produced a voltage 
difference across the middle layer. Leads of the EP thermoeleisent were spot- 
welded to tabs on the top and bottom layers. 

The manufacturer listed the uncertainty In the calibration of the meters 
as ±5%, with a linear relationship between heat flow and measured voltage 
output. As used In this application, though, they were not interchangeable. 
That Is, If the plate was positioned such that each meter was exposed to the 
same heat flow, the output varied between meters by as much as 20%. This 
variability was not due to an uncertainty of an individual meter calibration 
constant, but rather to differences In the calibration constant between 
meters. Ve found Individual meters to be linear and repeatable to better thar, 
±1%, over a range of heat fluxes from 0 to 44 V/cm^ , from calibrating the 
meters In a simple guarded hot plate apparatus. The calibration constant was 
also estimated with the meters attached to the calibration specimen. Using 
the heat absorbed In the cooling water as an estimate of the heat flow through 
the meter, the linearity for a single meter was about ±2.5% for heat fluxes of 
18 to 52 W/cm^ , which was the same as the uncertainty In the measurement of 
heat absorbed In the cooling water. 

Ve believe differences In the calibration constants between meters are 
caused by: (1) variations In the thickness or thermal conductivity of the 
BKter; (2) differences In the absorptivity between the meters; (3) variations 
In the thermal contact between the meter and copper plate for different 
meters. Ve think that (3) Is the most likely; In the guarded hot plate 
apparatus we found that the agent bonding the meter to the cold sink affected 
the meter calibration. For exaiqtle, the calibration constant differed by 15% 
for a meter soft -soldered to the cold sink and for a meter with indium 
compressed between the meter and cold sink. 

As used In the furnace calihretlon, we obtained a relative calibration 
by assuming each meter was exposed to the same heat flow when located at the 
center of the furnace target. The voltage obtained at the center location was 
then used to normalize readings at locations off center. Based on the 
linearity and repeatability of the meters, the basic accuracy of the 
normalized reading was taken as ±1%. The uncertainty In the over-all heat 
flux distribution Is evaluated below. 
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A. 3. Methods of dat? analysis 

The heat flow neters measured the heat flow over a 2.32 cm^ area. For 
the specimen heat transfer experiments, heat flux (heat flow over an 
infinitesimal area) is required as a function of x position, y position, and 
the heater lamp voltage. It is more convenient to express this as 

q(x,y,v) = f,(x,y,v) Qj/A„ , (A.l) 

where q(x,y,v) s local heat flux based on the normal area; 

A„ - specimen normal area; 

V - heater voltage; 

fq(x,y,v) s heat flux distribution function. 

The quantity Qt/Aq is the average heat flux over the specimen. As re-stated, 
the problem now is to find f, . The steps in determining the distribution 
function are 

1) Measure the heat flow meter output (heat flow) over the target area 
(vs X and y) for several heater voltages. 

2) For each voltage, normalize the readings by the value of the meter 
output at the center of the target. 

3) Fit a curve to the normalized data vs. x and y. 

A) Using the fitted curve, transform the data to extract the local 
normalized heat flux. 

These steps are now explained in more detail. 

Thr meter heat flow is equal to the measured voltage times a calibration 
constant, K, or: 

jt(x.y.'v) “ K V(x,y,v), (A. 2) 

with jfc(x,y,v) = heat flow; 

V(x,y,v) a the measured meter voltage. 

B}' normalizing all the readings by the reading at the center of the furnace 
(x=0,y«0), we can define another function: 

J(x,y,v) * V(x,y,v)/V(0,0,v) * (x,y,v)/Jt (0,0,v) . (A. 3) 

This fvinctlon, through the linearity of the meter, is also the ratio of heat 
flow at x and y to the heat flow at the furnace center, tflth values J(x,y,v) 
at a few discrete points, typically a scan of x at constant y or a scan of y 
at constant x, we fit a curve to the data using a least squares regression 
scheme to develop a polynomial function for j(x,y,v). The data show that the 
variation in j with y is less than the uncertainty in the heat flow meter, so 
we neglect it and determine only a dependence of J on x. For simplicity in 
the notation, the v argument is dropped also, with the understanding that the 
operations which follow occu'* at a constant voltage. 
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tfe now define « function f(x) which ie Che local heac flux noraalized to 
he heat flux at x=0. 

f(x) ^ q(x)/q(0). (A-4) 

The function j ia Chua the integration of f over the awCer area, or 

J(x)-2w = P’f(f)df . (A. 5) 

Js-m 

where w ia the half-width of the awter. If f ia conatant or linear over the 
■eter area, eq (A. 5} will ahow chat J f. The relationahip between f and f,, 

la found by coad>ining (A.l) and (A.4); 

f(x) q(0) = q(x) ^ f,(x) Qt/A^. (A. 6) 

By integrating (A. 6) over the total target area and equating to , it can be 
ahown that; 


f,(x) = 


f(x) 



+Xq 

(l/(2xo)) , f(f)df 


(A. 7) 


where Xq la the half width of the target. Hence if we can find f(x) through 
eq (A. 3), eq (A. 7) can be used to find the diatribution function f. and eq 
(A.l) ia uaed to find q(x). 


To find f(x) , we define a function g such that 
J(x)-2w = g(x+w) - g(x-w) = r d(g(f)). 

J* - W 

Proa eqa (A. 5) and (A. 8} we see that 


(A. 8) 


f(x) d(g(x)]/dx. (A. 9) 

Considering for the aoiaent x>0, froa eq (A. 8), we can substitute *=x4w. and 
then 


g(z) * g(z-2w) + 2wj(*-w). (A. 10) 

ihe function g at z-2w can be evaluated by recursion of eq (A. 10), or 

g(z) * 2wJ(z-w) + 2wJ(z-3w) ♦ g(z-4w). (A. 11) 

And for an expansion of n teras, 

n 

f(z) « g'(z) * 2w- 53 ♦ f(z-2nw), (A. 12) 

1-1 
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where cq (A. 9) has been used where appropriate. The nuiiber of terns n is 
chosen such that ~w<z-2nw<w; that is, f(z-2nw) is evaluated around the center 
of the furnace target. r*se measurements show that j(x) is nearly linear near 
the target center, and f can be taken equal to J with little error. Or, 
replacing z with x. 


n 

f(x) = 2w 2^ (j' (x-(2i-l)w] I + j(x-2nw). (A. 53) 

1=1 


To evaluate f at a location x, one suaa up J' evaluated at x-w, x-3w, ... 
and adds to j near x=0. It Is easily shown that for x<0, 

n 

f(x) = -2w- £ (J’(x+(21-l)w)) + J(x+2nw), (A,U) 

1=1 


and n is chosen so that -w<x+2n»Kw. 

In suonary, j(x) is calculated fron the neter voltages for the aeasured 
points fro« eq (A. 3) and fitted to a pTlynonlal . The function f(x) Is 
calculated fron eqs (A. 13) and (A. 14), and also fitted to a polynonla) 

Finally, f. Is calculated fron eq (A. 7). The najor advantages for <ising this 
technique (as opposed to using J(x) directly) to detemlne the calibration 
function are two- fold. First, If the actual heat flux varies with x In a non- 
linear nanner, the neasured J Is not equal to f. This occurs near the end 
walls In the furnace, as the results rf the next section will show. Second, 
the nethod allows for extrapolating to the furnace wall (x=±Xg), even though 
the heat flow neter, when aligned with the wall, gives a heat flow centered at 
distance w fron the wall. This Is because J ' (x) represents the change Ir heat 
flux between (x+w) and (x-w), and j at a neter half-width fron the wall gives 
f at the wall. 

A. 4. Results of the calibration 

Ve calibrated the furnace at a range of heater lanp voltages fron 201 of 
full voltage to 50% of full voltage. Very roughly, the heat flux on the 
target la linear with heater voltage; 30% of full voltage corresponds to about 
30 V/cn* of heat flux on the target. At a single calibration, we held the 
heater voltage constant and neasured heat flow through the neters over a 
representative portion of the furnace target. Voltages below 20% wc.^ ..ot 
calibrated as they corresponded to heat fluxes lower than necessary for the 
heat transfer experlnents, and voltages above 30% were not calibrated to 
avoid overheating the heat flow neters. The calibration testa are sunnarized 
in table A.l. We repeated the tests at voltages of 20%. 33%. and 30% to 
aaseaa Che reproducibility of the data. At the 33% voltage level, we 
conducted calibrationa at several levels of air cooling of the heater to 
detemlne whether Che cooling rate had an affect on the calibration. 

Ve present results in terns of the neasured nomallzed heat flow through 
Che heat flow neters (J of Che previous section) plotted as a function of the 





position variables x and y. This indicates the trends of the data. C’^rve 
fits of the heat flow data are plotted where appropriate, and the calibration 
equations ( f ^ function) for the dimensionless heat flux are presented. 


Table A.l. Summary of furnace calibration experiments 


Test 

Nunber 

Date 

Heater 
Voltage .. _ 

Type of 
Scar» 

CooBBents 

1 

6/13/88 

35.3% 

x-scan 

Scan 

at y-0 and both 





side 

walls . 



35,3% 

y-scan 

Scan 

at x=0 and both 





end walls 

2 

6/22/88 

35.3% 

x-scan 

y=C; 

3 cooling air 





flow 

sett ings . 




y-scan 

x=0. 


3 

6/23/88 

19.9% 

x-sca;i 

y=0. 





y-scan 

x=0. 


4 

6/23/88 

50.8% 

x-scan 

y=0. 





y-scan 

x=0 


5 

8/16/88 

35,4% 

X'sean 

y-0; 

2 cooling air 





flow 

settings . 

6 

8/18/88 

20.4% 

x-scan 

y=0. 


7 

8/18/88 

35,2% 

X - scan 

y=0. 


8 

8/18/88 

50.0% 

x-scan 

y=0. 


9 

9/06/88 

20.2% 

x-scan 

y=0. 


10 

9/06/88 

24.5% 

x-scan 

v=0. 


11 

9/06/88 

29.3% 

x-scan 

y=0 


12 

9/06/88 

35.2% 

x-scan 

y=0 



Traverses in the y (short) direction of the furitace indicated uniform 
heat flow at all voltages tested. This can be seen in figure A 3, which 
shows the normalized meter heat flow (J function) from testa 2 and 3, voltages 
33% and 20% respectively. Tie scan was done near x=0. the midpoint of the 
furnace. Data from meters 1 and 2 are shown; the data overlap near y^O. 

There is about a 3% drop off in the meter output near the side-walls (±2.3% of 
the mean), but since this is of the sane order as the meter calibration 
uncertainty, we neglect this effect in developing the calibration functions. 

)n test 1, the furnace was scanned in the y direction clos'i to both end walls 
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(large x and small x) . The results indicated the same uniformity in the heat 
flux for the y direction as in figure A. 3. Hence, we found no significant 
variation in heat flux in the y direction, regardless of heater voltage or x 
position in the furnace. 

Scans in the x (long) direction showed variations with position and in 
some cases the heater voltage. Figure A. 4 shows meter heat flow data for 
tests 2, 4, and 11 using heat flow meter #3, at voltages of 35%. 51%, and 29%. 
This distribution was unsynmetric; the heat flow gradually increased from the 
negative x end to the positive x end, with a drop-off near both end walls. We 
measured this asymmetry for all voltages above 29%, which we believe is due to 
the air cooling characteristics of the radiant heater. In the heater, cooling 
air flows around the lamp elements to prevent overheating; the cooling air 
flow direction is from the negative x end to the positive x end. In flowing 
along the lamps, the air heats up significantly. Hence the cooler lamps at 
the air inlet (negative x) produce lower heat flux that the hotter lamps at 
the air outlet. Reversing the position of the heater relative to the furnace 
produced reverse asymmetry in the heat flow distribution. 

A curve fitted to the data from tests 1, 2, and 4 is also shown in figure 
A. 4, and the standard deviation of the residuals between the measured heat 
flow and the curve value is 0.5%. Although the heat flow varies with x 
position, we found no significant change in the distribution above 29% 
voltage. When data from later tes^s at voltages above 29% (tests 7, 8, 11, 
and 12) was compared to the curve-fit from figure A. 4, the standard deviations 
to the curve for these later tests ranged from 0.9% for test 11 to 1.3% for 
test 7, indicating the repeatability of the furnace. In test 2 we varied the 
cooling air flow rate such that the temperature of the air exiting the heater 
changed from to 250*^0 (lower cooling air flow corresponding to higher 

air temperature; . This had no significant effect on the heat flow 
distribution. 

For heater voltages of 25% and lower, the furnace produced a symmetric 
heat flow distribution in the x direction. Figure A. 5 shows scans in the x 
direction from tests 6, 9, and 10 for heater voltages of 20% and 25%. Also 
shown is the curve fitted to the meter heat flow data (j function) , which had 
a standard deviation of 0.3% from the measured data. The heat flow drops to 
about 85% of the maximum at the end walls. The variation in the y direction 
was again no gre ter than that shown in figure A. 3. We believe that at the 
low voltages, the heater lamps remain cool along their entire length, which 
results in the symmetric distribution. 

The measured heat flow data for tests 1, 2, and 4 were analyzed to give 
the furnace distribution function (f^) for the unsymmetric profile (voltages 
greater than 29%), and heat flow data from tests 6, 9, and 10 were analyzed to 
give the distribution function for the symmetric profile (voltages less t tan 
29%). The sixth-order polynomial curve-fit and the region validity are 
summarized in table A. 2. The actual heat flux at any x an., y for a heater 
voltage v is then given by eq (A.l). 
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Table A. 2. Furnace heat flux distribution function, fq 


Type of Profile 


Region of Validity Coefficients of 

PoWnoirial, 


Synunetric Lamp Voltage < 29% Aq = 



Unsymmetric Lamp Voltage > 29% Aq = 

Ai = 

A2 = 
A3 = 
A4 = 
A5 = 
Ag = 

* where f^ = Aq + A^ x + A 2 -x^ + A 3 x^ + A^ *x^ A 5 
and X is in cm. 


1.053 

205 


6.291 

236 

E-5 

-1.770 

100 

E-3 

-6.021 

114 

E-5 

4.450 

009 

E-5 

1.028 

527 

E-6 

-1.750 

378 

E-6 

1.035 

180 


2.761 

752 

E-3 

-1.480 

886 

E-3 

4.125 

198 

E-4 

7.524 

066 

E-5 

-8.586 

921 

E-6 

-2.047 

184 

E-6 

+ Aj 

i-x® 



A. 5. Uncertainty in the furnace heat flux calibration 

The uncertainty in the heat flux distribution is estimated from the 
uncertainty in the heat flow meter measurement, a, the uncertainty in the 
calibration function when compared to the measured heat flow data (3a), and 
from the reproducibility of the furnace calibration (3a). The square root of 
the sura of the squares of these individual uncertainties results in an 
uncertainty of about ±4%. 
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1.3 mm 












Figure 3. Speclaen furnace. Inlet and outlet gaa teBperatures are 

measured at points 0 and 1; upstreaa and downstream pressures 
are measured at points A and B. 
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Thermocouple 



Glass Siesvs. 

0.51 
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NicrosI 
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0.25 


NOTE 

All dimensions in mni 


Figure 4. Technique for aountlng cheraocouples on pin- fin speclnens. 
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Figure A. 2. Lay-out of heat flow Mtera on calibration specinen 







Normalized Meter Heat Flow 
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Figure A. 3. Traverse of normalized heat flow along y direction at x»0 
Shorn are data from tests 2 and 3. heater voltages 35.3% 
and 19.9%. **" is meter 1, and *o” is meter 2. 
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Figure A. 4, Traverse of norsMlized heat flow aloAg x direction at y=0 
for heater voltages greater than 29%. Shown are data froa 
test 2 (35.3% voltage, •*"), test 4 (50 8% voltage, "o"), 
and test 11 (29 3% voltage. The saooth curve is the 

curve fitted to the noraalized heat flow data (J function), 
which has a standard deviation of 0.5%. 
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Normalized Meter Heat Flow 





X Position (cm) 


Figure A.S. Traverse of norflulized heat flow along x direction at y=0 
for heater voltages less chan 29%. Shotm are data froa 
teat 6 <20.4% voltage, test 9 (20.2% voltage, "o"), 

and test 10 (24.5% voltage, "dk”)* The saooth curve is the 
curve fitted to the noraalized heat flow data (J function) , 
which has a standard deviwtion of 0.3%. 
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